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We present an investigation of a gauge invariant smearing technique that allows the construction 
of smearing functions with arbitrary radial behaviour, by foresaking the space filling nature of 
traditional smearing techniques. This is applied to both heavy-heavy, heavy-light, and light-light 
systems with one particular choice of radial "wavefunction" - the hydrogenic solutions - and we find 
ON ' good stability for both fitted masses and amplitudes of the radially excited states. The dependence 

of the amplitudes on the smearing radius is demonstrated to be well understood, while near optimal 
smearing radii may be found with extremely low statistics using a property of the smeared-local 
correlator. The smearing technique is inexpensive since it is non-iterative, achieves a good signal 
S*H , to noise ratio, and can be altered to use wavefunctions from, say, potential models or the Bethe- 

■ Salpether equations in future simulations. 

<n : 
<n . 

I. INTRODUCTION 

mi 

Smearing is a technique that has been used for some time in Lattice QCD simulations |lp|3|] in order to improve the 
quality of the results obtained for a given computational expense by reducing the extent of excited state contamination. 
The choice of operator used to represent a state in any given simulation is essentially free within the constraint that 
it must have the appropriate quantum numbers of the desired state, which (usually) guarantees non-zero overlap 
, with the groundstate of those quantum numbers. In Euclidean space, all such operators asymptotically behave as the 
lowest lying state with quantum numbers matching those of the operator J]. We consider the two point correlation 
function of a set of operators Ok with some definite J PC in a Euclidean lattice of temporal extent T. Inserting a 
'"" ■ complete set of physical intermediate states of the same quantum numbers we find the same asymptotic behaviour 
Q_i| for all operators 

|: (o k ( x) ot m = E <o|0fcl ;l n|0i|o) ^-co S h^ft-|) a) 

= ^2C kn C* kn ^^co S hE n (t-^) (2) 

c3 ; n ™ 

-E„f T 

This somewhat egalitarian situation is the principal reason why lattice calculations have traditionally had difficulty 
extending their scope to radially excited states. However, some choices of operators turn out to be more equal than 
others, depending on the relative values of the C n k- For ground state phenomenology we would like Cok C n k Vn > 
1. To this end it is common to construct extended operators which project more predominantly onto the groundstate 
of the desired system; it is this extension procedure that has been dubbed smearing. Typically we require that the 
smearing procedure be parity and charge positive, and possess at least cubic symmetry 

There are two different schools of thought on how to construct smeared operators - gauge invariant smearing and 
gauge fixed smearing. 



* Present Address: Department of Physics and Astronomy, University of Glasgow, Glasgow 

1 Ignoring questions such as scalar glueball mixing in quenched simulations, and as we shall see special cases where the overlap 
is fine tuned to be near zero. 



1 



A. Gauge Invariant Smearing 



Explicit gauge invariance is one of the most attractive features of non-perturbative simulation techniques, and 
remains to this day a necessary (but sadly insufficient) condition for correct operation of simulation code. Further 
the gauge fixing process may be multiply defined due to Gribov ambiguities of some gauge fixing conditions j|,[3| , and 
one has to assume that the ambiguity does not cause any bias in the quantities being measured. For these reason 
many people are unwilling to use gauge fixing, and proceed with one of a number of gauge invariant techniques. 

The most commonly used gauge invariant techniques are Wuppertal H , Jacobi [0 , and fuzzed smearing || . Wup- 
pertal smearing covariantly smears the quark fields with a spatial distribution which is the solution of the scalar Klein 
Gordon equation on a given gauge configuration, while Jacobi smearing is a numerically efficient way to approximate 
this. In order to preserve gauge invariance both the first two methods iteratively apply a gauge invariant nearest 
neighbour operator to smear out a quark field, with limiting behaviour to weight sites in a "bell shaped" distribution. 
Given the extended nature of the physical states, and the wavefunction picture in the non-relativistic limit one would 
expect smeared operators to more accurately project onto the ground state, and indeed this is borne out by simula- 
tions. However these gauge invariant techniques give little control over the precise functional form of the smearing 
beyond a basic radius parameter. 

The cost of iteratively applying, say, the Jacobi operator N times on a propagator is somewhat prohibitive. Lacock 
and Michael j|] suggested the use of a covariantly transported "cross" using low noise fuzzed links. The fuzzing 
prescription iterates 

U'j(x) = c U 3 (x) + [U? tapl >, x +j) + U st r l °(x, x + j)] , (4) 



U j (x)=VU' j (x), 



(5) 



where V is a projection onto SU(3) (not closed under addition) by Cabibbo Marinari maximisation of Tr[[/|[/'] with 



six hits on each SU(2) subgroup. Typically c is set to 2.0, and there are 5 iterations performed. These fuzzed links 
are used to transport the quark field by some number of sites N along each of the spatial axes. 



i— x,y,z 



11 Ui{x + (n-l)i) 



71=1 



x ip{x + ri) + 



l[ul(x-rri) 



n=l 



x ip(x — ri) 



(6) 



The fuzzing radius, r, is chosen to minimise the contamination of the ground state operator. 



B. Gauge Fixed Smearing 

Gauge fixed smearing (either Coulomb or Landau gauge) allows complete freedom for the functional form of the 
smearing; in a fixed gauge we are freed from the need to iterate local covariant transport to space fill, and convolve 
quark fields with arbitrary smearing functions using fast Fourier transforms. In particular, the use of smearing 
functions with nodes has allowed reliable extraction of excited state energies from NRQCD through the use of multi- 
exponential fits [|| . Neither the cost of Fourier transforming the data set multiple times to implement this smearing 
on a massively parallel machine nor the loss of gauge invariance are desirable features, and there is clearly a need for 
a gauge invariant smearing technique with free radial functional form. 

In order to proceed, we consider the non-relativistic interpretation of qq systems in terms of spatial wavefunctions. 
This is certainly appropriate to heavy-heavy mesons, however, we can apply this technique in heavy-light and light- 
light systems without needing to worry about its relevance. 



II. SMEARING IN A NON-RELATIVISTIC POTENTIAL 



We consider smearing in the context of a non-relativistic potential model with spherical symmetry. The solu- 
tions may be written as ipnim(r, 0, (ft) — Rni(r)Pi m (6 ', 4>) where Pi m (9,4>) are the usual spherical harmonics. The 
orthogonality relation of the wavefunctions is then 
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dV^ ko ip n ip = 5 k iS op / R* ml (r)R n i(r)r 2 dr 



(7) 



r=0 



R ml {r)R nl {r)r i dr = 8 mn (8) 

We shall ignore the / and m indices from here, and consider spherically symmetric functions. Suppose we construct 
non-local meson operators 

O n (x) = J d 3 yi>(y)S n (y,x)ij;(x) (9) 

where the smearing function S n (y, x) — S n (y — x) is thought of as a spatial wavefunction defined on lattice sites 

S n { Vl Q) = S n (y)Y,S z (y - ix - jy - kz) (10) 

Note we have dropped the colour indices since in our wavefunction approximation the gauge degrees of freedom are 
represented by a colour singlet potential. 

The operator O will have an overlap with each physical state \ip m ) of the relevant J PC 

C nm = J ip* m {x)S n {x)dV = ]T ip* n (ix+jy + kz)S n {ix + jy + kz) (11) 

i,j,k 



Now if S n (x) is chosen to well approximate the state ?p n (x) and the lattice is sufficiently fine that the sum ([11]) is 
close to the integral 



Cnm = J r m S n {x)dV ~ / r m (x)^n(x)dV = S nm (12) 

Thus with a well chosen set smearing functions S n we can write 

Cnjn — &mn ~\~ Tj nra (1^) 

where r\ nm is a small contamination to the signal, giving a set of correlation functions corresponding to each of the 
physical radial excitations. In fact rj nm ~ Vra < n guarantees the lowest state significantly contributing to the 
correlation function of the nth radially excited operator is in fact the nth state, and the effective mass will not decay 
to the ground state of the J PC on a finite lattice. 

If we choose a non-space-filling smearing function of the form 

N 

S n (x,0) = ^r 2 0„(r) ]T (S^x-r^ + S^x + rfi)) (14) 

r— jl—x,y,z 

where <pn is some arbitrary modulating function that is chosen to approximate the true solution ip n , and the factor 
r 2 compensates for the non-space-filling nature of the smearing, then the overlap is 



Cnm = J r m (x)S n (x,0)dV (15) 
N 

= ]Tr 2 </>„(r)V4(0 (16) 



r=0 



r 2 r m 4>ndr (17) 
S nm - (18) 



r=0 
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Thus if we choose the 
need for space filling, 
colour: 



<j> m correctly (i.e. to be ip n ) then we establish the orthogonality relation required without the 
The new smearing technique proposed is a generalisation of the above form to incorporate 



N 



tp(x) -+^( r +i) 2 0„(r) 



Yl Ui(x+(n-l)i) 



n 



(19) 



where the links are fuzzed links, c.f. Section IA. We use the factor (r + ^) 2 as opposed to r 2 so that there is a 
non-zero contribution from the local current, since this was observed to improve the statistical noise. The procedure 
is very similar in cost to fuzzing with a radius N, since all the terms can be formed as part of the process of fuzzing 
to the largest radius. The technique therefore eliminates the two main problems with implementing gauge-invariant 
multi-state smearing, namely that gauge invariance requires expensive iterative space filling techniques, and that it 
restricts the functional form. Here we have a low-cost gauge invariant technique with an arbitrary functional form, 
allowing the insertion of both ground and radially excited wavefunctions. 

We shall consider the case when the physical solutions are approximately hydrogenic wavefunctions. 



III. OPTIMISING THE SMEARING 



We consider the correlation function with a smeared source denoted R and a local sink denoted L, and restrict the 
argument to one contaminating excited state, and ignore the periodicity of the lattice 



C t = Ae~ Elt + Be' E2t 



(20) 



where A oc CliC r1 and B oc CliC* R2 . The effective mass of this correlator is loosely the negative of the derivative of 
the log: 



Mi. 



dt 8 ' * 1 + B e -(E 2 - El )t ■ 



(21) 



By inspection it can be seen that there are two distinct cases for the approach to the plateau: 

• 3 > approach from above 

• ^ < approach from below, vertical asymptote at t = 

Suppose the physical wave functions are well described by the hydrogenic wavefunctions, ip n , with some Bohr radius, 
n. We wish to approximate these solutions with a trial wavefunction, which we choose to be the hydrogenic solutions 
4> m with Bohr radius T2- Consider the overlap of the trial excited state with the physical groundstate, 



Cv 



r ifjQtfiidr 



oc / r (1 )e 2r 2 e r i dr 

2r 2 



r(3) - — r(4) 



2c d - 

T2 



(22) 
(23) 
(24) 
(25) 



where c : 



2r±r 2 



S = Ti — 7*i, and the last step is made for small S. Thus the overlap of the excited state smearing 



with the physical state is proportional to the error in the Bohr radius. The overlap is positive for smearing radii that 
are too large, and negative for radii that are too small, and so the effective mass of the smeared - local correlation 
function will have a vertical asymptote for radii that are too small, and a smooth form for radii that are too large. 
This dramatic change of behaviour may be used to tune the Bohr radius on very small statistical samples. 

In fact it should always be possible to orthogonalise the first radially excited smearing with respect to the physical 
groundstate by tuning the Bohr radius parameter even when the physical states are not hydrogenic. One choice 
for the radius will not, however, simultaneously minimise the overlaps C±o and Cqi- The better the ansatz for the 
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smearing functions, the closer together we expect the optimal radii for the groundstate and excited state correlators 
to lie. Where the ansatz is better we also expect there to be less contamination from the higher excited states, since 
the physical wavefunctions are all mutually orthogonal. 

In order to demonstrate the applicability of this smearing technique, a feasibility study was carried out on a small 
sample, with a large number of smearing radii, in the light-light, heavy-light and heavy-heavy sectors. 



IV. FEASIBILITY STUDY 



We use the notation R n [ to denote smearing with the corresponding hydrogenic wavefunction with some Bohr 
radius. Fourteen quenched configurations of a 24 3 x 48 lattice at (3 = 6.2 were used with non-perturbatively 0(a) 
improved quarks (i.e. C sw = 1.61377). Kappa values 0.13460 and 0.12300 were used, corresponding to quarks near 
the strange and charm masses respectively. The lighter quark propagators had local sources, and both local and R\o 
sources were generated for the heavy propagator with Bohr radii of 2.4, 2.5, 2.6, 2.75, 2.8, 2.9, 3.0, 3.1, 3.25, 4.0, and 
5.5 in lattice units. This allowed the smearing combinations in Table | to be generated for each Bohr radius. The 
correlation functions were analysed using uncorrelated fits, because the small statistical sample caused the correlation 
matrix to be too noisy for reliable inversion in correlated fits. As such the estimates of % 2 are probably unreliable. 

Of primary interest is the simultaneous two exponential fit to the R\q and R20 source smeared - local sink correlation 
functions. We show the fitted effective masses for various Bohr radius for the heavy-heavy, heavy-light and light-light 
sectors in Figures |l| through |[ 

The effective mass for the radially excited smearing demonstrates the expected transition in behaviour, and appears 
to be forming a plateau for ro = 2.6(3.0,3.1) in the heavy-heavy (heavy-light, light-light) sector. While the formation 
of an excited state plateau is clearly strongly dependent on the selection of ro to within about 0.1 (this is possible 
to do using only one configuration by looking for vertical asymptote in the effective mass plot), we do not wish the 
fitted values to be dependent on the smearing. To investigate this we study the stability of the fitted masses with the 
input parameters to the fitting procedure, namely the Bohr radius, and the timeslice range of the fit. 



A. Stability Analysis of Fitted Masses 



The fitted values for the ground state and excited state masses are given in Table || to Table |v|, and are found to 
be remarkably stable in the fit range typically agreeing within statistical error for all values of tmin > 3. 

The level of statistical error on the excited state mass is also particularly pleasing. The error on the excited 
light-light state is typically between only 5 and 10 times larger than that on the ground state. 

In order to study the stability of the fits with respect to the chosen Bohr ra dius , I tabulate the fi tted mass values 
corresponding to t min — 4 for all the simulated Bohr radii in Table VI, Table VII, and Table VIII. They are stable 



within the statistical error over almost the entire range of radii simulated. 



B. Selecting the Bohr Radius 

From inspection of the effective mass plots presented, we can see that the optimal values of the Bohr radius 
parameter (at least for the R20 correlator) were about 2.6 for the heavy-heavy state, 3.0 for the heavy-light and 3.1 
for the light-light. Using the discussion in Section [Ol] we plot the fitted amplitudes of the local-smeared correlation 
functions versus Bohr radius, looking for the zero of the overlap of the R20 correlator with the groundstate. 

Figure [7, Figure |^ and Figure || present the dependence of the four fitted amplitudes on the Bohr radius used in 
the heavy-heavy, heavy-light and light-light systems respectively. The linear fits are a good approximation, producing 



the optimal values for the smearing radii given in Table IX 



While the effective mass plots presented were on fourteen configurations, it is clear that following the pole on the plot 
gave very accu rate optimisation of the Bohr radius, as demonstrated by later fits to the amplitudes. It was anticipated 



in Section III that, since our smearing function basis is not the physical one, the radius for minimal contamination 
of the groundstate by the excited state would differ slightly from the radii which gave minimal contamination of the 
excited state. This is not apparent from the data presented in Figure ^, Figure |^ and Figure ^| where there is no 
minimum of the Rio - excited state amplitude. We must bear in mind however that the contamination of the R\q 
correlator is measured at early times, whereas the contamination of the R20 correlator is found from its asymptotic 
time dependence, leaving the determination of the optimal Bohr radius from the Rio correlator much more subject to 
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contamination from higher excited states. For this reason, it is dangerous to look only at the local-smeared amplitude, 
since cancelling contributions from the higher states can cause a seemingly flat effective mass. 

Instead we study the smeared-smeared correlators, which have positive definite contaminating contributions, and 
minimise the amplitude for the contamination of the Rio - R\q correlator. These correlation functions were only 
generated for the heavy-heavy combination, and Figure ^ presents a sample fit to the pseudoscalar data. Here the 
contamination below timeslice 5 is apparent, as well as a weaker drift to the plateau beyond this point. It is this 
ultimate drift we must study. Figure [ll] plots the fitted contaminating amplitude for this correlator versus radius; we 
expect from Section HI a quadratic dependence on the smeared-smeared correlator. We obtain a minimising radius 
of ro — 2.1(4) with a quadratic fit. 



C. Use of Doubly Smeared Correlation Functions 



In this section I demonstrate that the signal obtained from doubly smeared correlators is useful. The double 
smearing is created by constructing the meson propagator from two propagators which are each source smeared, 
while, for the purposes of this demonstrate, we use local sinks. The double exponential fit to the double Rio and 
single i?20 smeared correlation functions are presented in Figure [l^ and Figure [l3] for the heavy-heavy and light-light 
systems respectively. 

Clearly this smearing combination, while not tailored for any particular state constitutes a valid signal that can be 
included in multi-channel multi-exponential fits, which is always of use, especially for the p-states. 



V. COMPARISON WITH FUZZING 



In order to fairly compare the level of statistical error with traditional techniques, I present double exponential 
fits to the i?io sink and the fuzzed sink, both with the optimal radius in Figure Q Figure [l5] and Figure [l(] for the 
heavy-heavy, heavy-light and light-light systems respectively, on the same set of local source propagators. 

It can be seen that the statistical error from the R\q correlator is much reduced with respect to the fuzzed data. 
Admittedly the fuzzed data does appear to plateau earlier, with a caveat however; the noise on the fuzzed data 
is such that there is quite possibly a "double approach" to the plateau, due to contributions from more than one 
excited state, as seen in the large statistics fuzzed data at (3 = 6.2 in Figure 17, so that with higher statistics the 
plateau may well not survive. It has also been found that the fuzzed-fuzzed combination tends to plateau after the 
fuzzed-local, indicating opposing contributions from excited states, and meaning that the plateau at early times is 
a misleading balance between contaminations of opposite sign (at least until after the fuzzed-fuzzed correlator has 
reached a plateau), c.f. Figure 17. 



VI. CONCLUSIONS 



This smearing method provides, for the first time, a gauge invariant technique for inserting arbitrary radial wave- 
functions in a smeared operator. So far results have only been presented using hydrogenic wavefunctions, though 
in principle any form could be used. Excellent stability of fits to both the ground and radially excited states was 
demonstrated, and the level of statistical error was significantly smaller than traditional techiques on the same en- 
semble. Tuning the radius of the wavefunctions allowed a clear plateau for the radially excited state to be isolated. 
An understanding of the dependence of the fitted amplitudes on the radius of the smearing was obtained, which can 
be used to determine more accurately the optimal radius. 

The feasibility study did not include a sufficient number of quark masses to extrapolate to the physical spectrum, 
however a rudimentary comparison of the 25* — IS* splittings with experiment is made in Table [x|, and can be seen to 
be plausible. The main source of error, particularly in the light-light system, is expected to be finite volume effects, 
since the Bohr radius required suggests the extent of the wavefunction is significantly bigger than the lattice used. 

Since the technique allows a free choice for the wavefunction it would be interesting to try to find a better basis 
than the hydrogenic wavefunctions. Harmonic oscillator wavefunctions have been tried and found to be inferior to 
hydrogenic. Other possible options are numerically solving the Cornell potential and inserting the solution as the 
smearing function, as has been performed by the SES AM collaboration in the gauge fixed context , or alternatively 
measuring the wavefunction within a lattice calculation and re-inserting this. It would be interesting to compute the 
i?2o smeared correlators on a large volume for light quarks since we expect finite volume effects on the current volume. 
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FIG. 1. Heavy heavy pseudoscalar fits for ro = 2.4,2.5,2.6. The plots contain the effective masses in lattice units for 
correlators with R20 (circles) and Rio (fancy squares) source smearings with local sinks. The fits are simultaneous double 
exponential over the timeslice range 4-20. A clear transition between a pole in the R20 effective mass and a continuous effective 
mass is seen as ro is tuned across its optimal value. 
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FIG. 2. Heavy heavy pseudoscalar fits r — 2.75, 3.0, 5.5 The plots contain the effective masses in lattice units for correlators 
with R20 (circles) and Rio (fancy squares) source smearings with local sinks. The fits are simultaneous double exponential over 
the timeslice range 4-20. 
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FIG. 3. Heavy light pseudoscalar fits, ro = 2.6, 2.8, 2.9. The plots contain the effective masses in lattice units for correlators 
with R20 (circles) and Rio (fancy squares) source smearings with local sinks. The fits are simultaneous double exponential over 
the timeslice range 4-20. 
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FIG. 4. Heavy light pseudoscalar fits ro = 3.0, 3.1, 4.0. The plots contain the effective masses in lattice units for correlators 
with R20 (circles) and 7?io (fancy squares) source smearings with local sinks. The fits are simultaneous double exponential over 
the timeslice range 4-20. 
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FIG. 5. Light light pseudoscalar fits, ro = 2.8, 2.9, 3.0 The plots contain the effective masses in lattice units for correlators 
with R20 (circles) and 7?io (fancy squares) source smearings with local sinks. The fits are simultaneous double exponential over 
the timeslice range 4-18. 
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FIG. 6. Light light pseudoscalar fits ro = 3.1, 3.25, 4.0. The plots contain the effective masses in lattice units for correlators 
with R20 (circles) and Rio (fancy squares) source smearings with local sinks. The fits are simultaneous double exponential over 
the timeslice range 4-18. 
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Quarkonium Smearing Amplitudes 
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FIG. 7. Fitted amplitudes for quarkonia correlation functions using the timeslice range 4-20. The contaming amplitude in 
the R20 source smeared correlator clearly shows a zero at ro ~ 2.6 
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FIG. 8. Fitted amplitudes for heavy-light correlation functions using the timeslice range 4-20. The contaming amplitude in 
the R20 source smeared correlator clearly shows a zero at ro ~ 3.0 
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FIG. 9. Fitted amplitudes for light-light correlation functions using the timeslice range 4-20. The contaming amplitude in 
the R20 source smeared correlator clearly shows a zero at ro ~ 3.1 
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FIG. 10. Rio - Rio source and sink smeared quarkonium correlator (fancy squares), used in a double exponential fit with 
the R20 source, local sink correlator (circles) with Bohr radius r = 2.4, 3.0 
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- Ground State Contamin'ation - 
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FIG. 11. Contaminating amplitude obtained from the Rio - Rio source and sink smeared quarkonium correlator. This 
smearing radius minimising this amplitude lies in the region of 2.1(4). 





FIG. 12. Fit to the doubly Rio source smeared with local sink correlator (fancy squares) and the R20 source smeared with 
local sink correlator (circles) for the heavy-heavy meson. The quality of the plateau suggests that the use of double smearing 
is worthwhile for additional information on the groundstate. 
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FIG. 13. Fit to the doubly Rio source smeared with local sink correlator (fancy squares) and the R20 source smeared with 
local sink correlator (circles) for the light-light meson. 
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FIG. 14. Heavy-heavy fuzzed (r=6) and Rio comparison. Both correlation functions are sink smeared with local sources 
on the same configurations. The lower statistical scatter our smearing technique is thought to be in part due to including a 
non-zero contribution from the local current. 
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FIG. 15. Heavy-light fuzzed (r=6) and Rio comparison. Both correlation functions are sink smeared with local sources on 
the same configurations. 
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FIG. 16. Light-light fuzzed (r=8) and Rio comparison. Both correlation functions are sink smeared with local sources on 
the same configurations. 
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FIG. 17. Examples of approaches to fuzzed plateaux with high statistics. Opposing overlaps with the second and third 
excited states can cause non-monotonic approaches. This is shown in the vector plot, where the open circles represent the 
fuzzed source, local sink state. The pseudoscalar plot shows a false plateau in the fuzzed-local correlation function, since the 
fuzzed-fuzzed correlator forms a plateau later we know that the early points in fuzzed-local correlator plateau cannot contain 
pure information about the groundstate. 
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TABLE I. Smearing combinations 
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TABLE II. Heavy heavy double exponential fit stability; ro = 2.6, T ma x = 20 





v 2 Irlnf 
X l ao J 


v 2 
X 


Pom l /"I /"i 

JT be LI LIU 


Pcoi i r\ t~\ 
JT be LI LIU 


2 


0.66 


21. 


1.123(3) 


1.45(2) 


3 


0.25 


7.5 


1.125(4) 


1.42(1) 


4 


0.13 


3.7 


1.127(4) 


1.41(2) 


5 


0.13 


3.3 


1.128(4) 


1.40(3) 


6 


0.13 


3.2 


1.129(4) 


1.40(3) 


7 


0.14 


3.1 


1.130(5) 


1.40(4) 



TABLE III. Heavy heavy triple exponential fit stability; ro = 2.6, Tmax = 20 



tmin 


X 2 /dof 


x 2 


Pseudo (2300,2300) 


Pseudo' 


Pseudo" 


2 


0.081 


3.6 


1.128(4) 


1.39(3) 


2.13(4) 


3 


0.083 


3.5 


1.129(5) 


1.39(3) 


2.10(10) 


4 


0.086 


3.3 


1.130(6) 


1.39(4) 


1.9(2) 


5 


0.089 


3.2 


1.130(5) 


1.40(7) 


1.9(2) 


6 


0.095 


3.1 


1.131(6) 


1.39(8) 


1.8(2) 



TABLE IV. Heavy light double exponential fit stability; r = 3.0, T max = 20 



tmin 


X 2 /dof 


x 2 


Pseudo 


Pseudo' 


2 


0.44 


14. 


0.746(5) 


1.18(2) 


3 


0.15 


4.5 


0.745(6) 


1.13(2) 


4 


0.074 


2.1 


0.744(6) 


1.10(2) 


5 


0.076 


2.0 


0.744(7) 


1.09(3) 


6 


0.078 


1.9 


0.744(7) 


1.08(4) 


7 


0.083 


1.8 


0.744(7) 


1.07(6) 



TABLE V. Light light double exponential fit stability; r = 3.1, T ma x = 18 



tmin 


X 2 /do.f 


x 2 


Pseudo 


Pseudo' 




2 


0.18 


5.0 


0.280(7) 


0.86(4) 




3 


0.085 


2.2 


0.279(7) 


0.81(3) 




4 


0.079 


1.9 


0.278(7) 


0.79(4) 




5 


0.082 


1.8 


0.278(7) 


0.81(6) 
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TABLE VI. Heavy heavy double exponential fit stability with Bohr radius 



Bohr 


tmin 




x 2 


Pseudo 


Pseudo ' 




2.60 


4 


0.13 


3.7 


1.127(4) 


1.41(2) 




2.75 


4 


0.083 


2.3 


1.128(4) 


1.40(2) 




Z.oU 


A 

t± 


U.UDD 


1 Q 








2.90 


4 


0.045 


1.2 


1.128(4) 


1.39(2) 




3.00 


4 


0.033 


0.92 


1.128(4) 


1.39(2) 




3.10 


4 


0.026 


0.74 


1.128(4) 


1.39(1) 




3.25 


4 


0.021 


0.59 


1.128(4) 


1.39(1) 




4.00 


4 


0.016 


0.44 


1.128(3) 


1.39(1) 




5.50 


4 


0.016 


0.44 


1.129(3) 


1.40(6) 





TABLE VII. Heavy light double exponential fit stability with Bohr radius 



Bohr 


tmin 


x7*>/ 


x 2 


Pseudo 


Pseudo ' 


2.60 


4 


0.059 


1.6 


0.745(6) 


1.13(3) 


2.75 


4 


0.10 


2.8 


0.745(6) 


1.12(2) 


2.80 


4 


0.12 


3.2 


0.745(6) 


1.12(2) 


2.90 


4 


0.12 


3.3 


0.744(6) 


1.11(2) 


3.00 


4 


0.074 


2.1 


0.744(6) 


1.10(2) 


3.10 


4 


0.045 


1.2 


0.744(6) 


1.10(2) 


3.25 


4 


0.026 


0.72 


0.744(6) 


1.09(2) 


4.00 


4 


0.013 


0.35 


0.744(6) 


1.09(3) 


5.50 


4 


0.014 


0.38 


0.745(6) 


1.10(5) 



TABLE VIII. Light light double exponential fit stability with Bohr radius 



Bohr 


tmin 


X 2 /dof 


x 2 


Pseudo 


Pseudo ' 


2.60 


4 


0.0053 


0.13 


0.278(7) 


0.79(4) 


2.75 


4 


0.010 


0.24 


0.278(7) 


0.78(4) 


2.80 


4 


0.013 


0.32 


0.278(7) 


0.78(4) 


2.90 


4 


0.029 


0.70 


0.278(7) 


0.78(4) 


3.00 


4 


0.077 


1.9 


0.278(7) 


0.78(3) 


3.10 


4 


0.079 


1.9 


0.278(7) 


0.79(4) 


3.25 


4 


0.023 


0.55 


0.278(7) 


0.79(4) 


4.00 


4 


0.0013 


0.031 


0.277(6) 


0.77(7) 


5.50 


4 


0.00072 


0.017 


0.277(6) 


0.8(1) 



TABLE IX. Optimal Bohr radii 



System 


Fitted Value 


Choice 


2300-2300 


2.62(2) 


2.6 


2300-3460 


2.98(2) 


3.0 


3460-3460 


3.09(1) 


3.1 
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TABLE X. Comparison of 2S-1S splittings with experiment using string tension scale 



System 


Lattice (MeV) 


Experimental (MeV) 


Light Pseudoscalar 


1370(50) 


1170(100) (tt) 


Light Vector 


1250(100) 


700(25) (p) 


Heavy-Light Pseudoscalar 


970(60) 




Heavy-Heavy Pseudoscalar 


710(60) 


590 (J/ip) 
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